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Abstract 

We show that, when graphene is only subject to strain, the spin connection gauge field 
that arises plays no measurable role, but when intrinsic curvature is present and strain is 
small, spin connection dictates most the physics. We do so by showing that the Weyl field 
associated with strain is a pure gauge field and no constraint on the (2 + l)-dimensional 
spacetime appears. On the other hand, for constant intrinsic curvature that also gives 
a pure-gauge Weyl field, we find a classical manifestation of a quantum Weyl anomaly, 
descending from a constrained spacetime. We are in the position to do this because we 
find the equations that the conformal factor in (2 + l)-dimensions has to satisfy, that 
is a nontrivial generalization to (2 + l)-dimensions of the classic Liouville equation of 
differential geometry of surfaces. Finally, we comment on the peculiarities of the only 
gauge field that can describe strain, that is the well known pseudogauge field Ai ~ U 11 — U 22 
and A 2 ~ U 12 , and conclude by offering some scenarios of fundamental physics that this 
peculiar field could help to realize. 
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1. INTRODUCTION 


Graphene is a very promising table-top laboratory to indirectly probe some of the 
fundamental mysteries of Nature [1]. The low energy regime of its vr electrons is very 
well described by an effective theory that shares many of the features of a massless 
Dirac quantum held theory (QFT) in the presence of a background spacetime. 

In order for graphene to keep its promises, we need to have full control of what 
sort of helds are there and what they represent in a held theory language. In the vast 
literature on the gauge helds of graphene (see [2]-S], and [3] for a recent review), there 
are a variety of proposals, sometimes practically valuable for the applications to con¬ 
densed matter physics, but most of the time unsatisfactory for probing fundamental 
physics. The landscape of proposals ranges from SU{2) monopolelike gauge helds 
in the case of graphene membranes with intrinsic curvature (the inhated graphene 
buckyballs of [6]) to a concurrence of a spin-connection held and a U{1) held, in 
the case of purely strained graphene [7] (although sometimes non-Abelian helds are 
evoked in this case as well m)- Even in the simplest case, that is purely strained 
graphene, there is some confusion: does the spin connection arising from straining 
graphene give physical ehects or not? And, what is the interpretation of the U{1) 
held from a fundamental point of view? 

In this paper, we try to clarify part of this confusion, having in mind to come back 
to the full scenario of gauge helds in graphene in a forthcoming work. Our principal 
conclusions are the following: (1) When only strain is present, spin connection plays 
no role; only the U{1) pseudogauge held gives physical ehects. (2) When intrinsic 
curvature is present, spin connection plays a crucial role, to the point that, for small 
strain, it is spin connection that dictates most of the physics. In addition, we also 
clarify why some confusion might arise in relation to the non-Abelian structure (that 
is the Lorentz group) behind the Weyl (Abelian) pure gauge held [U HU] . 

To do that, we proceed as follows. We start oh with the fully relativistic approach, 
and proceed by seeing that the Weyl held associated to strain (that is what the 
spin connection reduces to in this case) is a pure-gauge held, that can never give 
a physical ehect, unless non-trivial topology is involved. Since we are in a pure 
strain conhguration (elastic deformations), and we only focus on local properties (our 
samples here have inhnite size), non-trivial topology of the Weyl held is excluded 
a priori. Later we show, for the hrst time explicitly, that this fully relativistic 
approach is equivalent to the standard (non-fully relativistic, static Hamiltonian) 
approach. We do so by showing that all the terms one obtains through the latter 
method are there also in the former. Having the equivalence, this closes the debate 
on the physical ehects of the spin-connection for pure strain. 

On the other hand, a crucial issue arises for the meaning of Weyl symmetry 
itself. In fact, the Weyl held is a pure gauge also when the membrane has nonzero 
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constant (hence non-trivial) intrinsic cnrvatnre, and we know that physical effects 
emerge from there HH. How can we reconcile this, on general gronnds, with the 
previons resnlt? We have solved the issne by first finding the eqnations that the 
conformal factor of the (2 + l)-dimensional theory has to satisfy, in the case of zero 
cnrvatnre (pnre strain) and in the case of constant cnrvatnre. This is interesting on 
its own right, as we have fonnd a nontrivial generalization to (2 + 1 ) dimensions of 
the classic Lionville eqnation of two-dimensional differential geometry of snrfaces, 
and of (1-|- l)-dimensional conformal field theories (CFTs). We then realize that the 
key difference between the flat and curved case, although they both give standard 
Weyl invariance, lie in the absence or presence, respectively, of constraints on the 
(2 -|- l)-dimensional spacetime. Therefore, in the curved case, the theory effectively 
lives in (1 -|- 1) dimensions where Weyl symmetry is anomalously broken. This is an 
interesting classical manifestation of a quantum Weyl anomaly. 

The paper flow is as follows: In Section we recall the results of previous work 
that we need here, and set the notation. In Sectionj^we find the connection between 
spatial and spatiotemporal conformal factors that allow us to conclude that the spin- 
connection cannot give physical effects for pure strain. We focus the last part of 
this section on the nonzero constant curvature case. We leave to Section 0] to show 
explicitly how the standard Hamiltonian and strain tensor language can be recovered 
from the fully relativistic formalism. After this top-down approach, in the Section]^ 
we use the tight-binding condensed matter Hamiltonian description, and comment 
on the peculiarities of the only gauge field that can describe strain. This latter road, 
that we call bottom-up, is actually necessary because, this kind of gauge field could 
not be guessed within standard QFT in curved spacetime. In the concluding Section 
we also offer some scenarios of fundamental physics that this peculiar gauge field 
could help to realize. 

2. SPIN-CONNECTION AND WEYL GAUGE FIELD FOR STRAINED 
GRAPHENE 

As is well known (see, e.g., m), the physics of the large wavelength vr electrons 
of graphene, at half-filling, can be efficiently encoded within the Dirac massless 
two-dimensional Hamiltonian. In our approach, we include time to make it fully 
relativistic, although with the speed of light c traded for the Fermi velocity vp (see, 
e.g., |9]); hence, we start from the (2 -|- l)-dimensional action 



( 1 ) 


where = (f, x, y) are the fiat spacetime coordinates, 7 “ are the usual Dirac ma¬ 
trices in three dimensions, and we expand around only one of the two Dirac points. 
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as we shall focus on strain. 

As shown in [H], the metric 


9^lu{q) 


1 0 0 

g 9ij i .^1 2/) 


( 2 ) 


can also describe strain; hence, we shall use the customary Dirac action in that 
curvilinear background. As also pointed out in [U], we recall that although the 
system is (2 + 1) dimensional, the Riemann tensor A,... = 0,1,2, has only 

one independent component, proportional to the Gaussian curvature, /C. Further¬ 
more, surfaces of zero or constant /C, make the metric ([^ flat or conformally flat, 
respectively, and both cases can be treated at once within a formalism that uses 


9i.u{Q) = 


( 3 ) 


where rjpy = diag(l, —1, —1), and the information about the metric being flat or 
not is encoded in the conformal factor S. The coordinates where the metric can be 
explicitly written in a conformally flat fashion. 


g^ = (T,X,F), 


( 4 ) 


are, in general, different from the original coordinates q^. 

Therefore, the natural candidate action to describe strained graphene is (for a 
while we shall use h = 1 = vp) 




where = ^u^^Jab, with Jab = |[7a,76], the Lorentz generators, is the inverse 
of the three dimensional vielbein e“ (the dreibein), is the spin connection, and 
being in (2 -|- 1) dimensions, we can write the spin connection in the one-index 
notation 



( 6 ) 


Here, a is the non-Abelian index of the local Lorentz group SO(2,1), and fi is the 
vector index on the spacetime manifold. In this formalism, the index a, a tangent 
spacetime index, plays the role of an internal/gauge index. 

The metric ([^ can always be written in more suitable coordinates = {t,x,y), 
where t is the same time for both coordinate systems, and {x, y) are the spatial 
isothermal coordinates of the surface 
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( 7 ) 



The single scalar function a identifies the surface/graphene membrane. Other 
isothermal coordinates can be found, say x, y, but the function identifying the surface 
is always the same^: 

a{x{x, y),y{x, y)) = a{x, y) . ( 8 ) 

3. CONNECTING SPACE AND SPACETIME CONFORMAL FACTORS 

When the surface described by a has zero or constant curvature, the two metrics 
([^ and Q both describe the same spacetime, although with different coordinates, 

qAQ) = ; ( 9 ) 

hence, their Gaussian curvature cannot differ. Due to the simple structure of the 
metric in the two coordinate frames, it is easy to compute the Gaussian curvature 

/C = -e-2'^[VV] (10) 

= e-2^[2nS + (4S)(a“S)] (11) 

where = 9| + 9?, and D = — dx — dy- Glearly, the two conformal factors are 

related, D((j): if we know Q^{q), we can write S(5;, y) = J2(T{x, y),X{x, y), Y{x, y)), 
and then knowing a{x,y), we obtain D((j). Nonetheless, we have the general equa¬ 
tions that S has to satisfy for the two cases, for /C = 0, 

□s = -i(a„s)(a“s), (12) 

which corresponds to a harmonic functions (i.e., solutions of = 0), and for 
/C = const ant 

□E = -^(a,S)(a“S) + i/Ce2^ (13) 

which corresponds to a Liouville functions (i.e., solutions of V^cx = —ICe^^). 


^ As an example of formula (|^, the metric of a sphere of radius r = 1 could be written as 
ds^ = {dx^ + dy^) with conformal factor a{x,y) = ln(l/coshy) being x S [0, 27r] and y G 

arcsin(l/ coshy), if y > 0 
TT — arcsin(l/coshy), if y < 0 
and X = X, a{x,y) = In sin y, with y S (0,7r), we recover the standard metric of the sphere 
ds^ = sin^ ydY + dy'^. 


(—oo,+oo). On the other hand, using the coordinates y = 
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Let us focus on the flat case which, since it corresponds to pure strained graphene, 
is the one of interest here. In this case, besides the obvious constant solution of (12), 
= C, we also have 


-‘flat 


= - -X^ -Y^) + C . 


(14) 


The constant C could be set to zero, but we shall keep it, to later compare with the 
curved case. For the conformal factor (12), the metric (|^ reads 


9t,AQ) = = 


^2C 


(T2 - X2 - F2)2 



(15) 


Hence, nothing constrains the norm of vectors just as for the Minkowski case 

„2C 


WQr = 9 ,uQ^q^ = 


(16) 


Vf.uQf^Q'' 

that means that truly describes a flat three-dimensional spacetime. In order to 
see explicitly that (16) is in fact the Minkowski spacetime, take the following change 
of coordinates: 


t = 


X = 


y = 


e^T 


J'2 _ _ y2 ’ 

'J'2 _ _ y'2 ’ 

e^Y 

J'2 _ Ji^2 _ y2 ■ 


( 17 ) 


( 18 ) 


On these coordinates the line element is ds^ = dt‘^ — dx'^ — dy^, showing that the 
singularities appearing in the metric (16) are simply coordinates singularities due 
to our choice of nonstandard coordinate. As a result of that, when we use T,fiat in 
the spin connection of ([^ 




( 19 ) 


we can safely use as a proper dreibein, because it indeed connects the tangent 
space with a flat manifold. Here and = claS, and we used the result 

that in three dimensions 7“Ja6 = 7b- 


2 To see whether one has a true or a coordinate singularity is, in general, not an easy task. On 
this, see, e.g., [T^ . 
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3.1. Zero curvature: no physical effects of strain through the spin- 
connection 


Therefore the action ([^ for the metric ([^ is 

.A = i y d^Q ^Y{da + ^a)i’ , 

and when the Dirac held is properly transformed 

xjj = , 

where refers to the Minkowskian hat spacetime [see 
simply 

A = i f , 


( 20 ) 


( 21 ) 


, the action (O 


IS 


( 22 ) 


which refers to the background metric 77 ^, 0 , which in turn is the unstrained situation. 
Strain is gone altogether. It has no physical effects. 


Another way to see this, is to stop at the action (20), i.e., before implementing 
the transformation of the spinor as in (21), and consider the gauge held Sq. This 


gauge held is itself a pure derivative, hence it cannot produce any physical effect 
through its held strength 


Fat = da^t - dt^a = (Saft ” SA)S = 0 , 


(23) 


as can be seen also by explicitly computing = V x T^fiat and = — VS^ + 


dx^fiat with Tjfiat in (14). The result is zero = 0 = 

Therefore, from here, we see that the very well-known pseudomagnetic held (and, 
for what matters, even a putative pseudoelectric held) induced by pure strain, cannot 
be accounted for by the spin-connection/Weyl pure-gauge held. 

Let us now comment on As seen, this is a pure gauge held associated to the 
local Weyl transformations (§ and d^ . Indeed, the Weyl held transforms as 

(24) 




The reason why we do not have here the full Weyl gauge held, IT^, but only its pure 
gauge part, is due to the local Weyl invariance of ([^; see [S]. Like any other Weyl 
held, Sq is an Abelian gauge held. Abelian gauge helds are those routinely used 
in the graphene literature, [3] and [7]. On the other hand, carries information 
on the non-Abelian structure of the local Lorentz group that is encoded in the 
spin connection. This information is in the tangent space index “a” of see (|^, 
(19) and discussion in between. Indeed, (local) Weyl transformations, in general. 
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rephrase spacetime scaling as an internal transformation. Non-Abelian gauge fields 
have also appeared in various discussions on the gauge held approach to strained 
graphene [8]. 

The above mentioned properties are common to the full Weyl gauge held 
hence, they hold also for theories that do not have local Weyl invariance. The 
extra property of is that it is SqS; i.e. the true degree of freedom is just one 
scalar, S, the one related to the two-dimensional spatial phonon, a, of the graphene 
membrane, as shown before. 

These facts are fully transparent in the coordinates in (|^. When the co¬ 
ordinates are diherent, the three aspects of this gauge held—(i) scalar nature, (ii) 
abelian held, and (hi) non-Abelian Lorentz structure—get mixed together, and they 
may appear, in the standard coordinates/frames used in graphene, as originating 
from diherent gauge helds, as we show in Section 


3.2. Nonzero curvature: the classical manifestation of the quantum Weyl 
anomaly 


Apparently, all seems clear: when E = Efiat, which should describe pure strain, 
no physical ehects can be described by the QFT in the curved spacetime approach. 
Nonetheless, we also saw that when the Gaussian curvature of the membrane is 
constant, a similar procedure could be applied; hence, this seems to lead to conclude 
that also in that case as well, there is no physical ehect. However, as we shall now 
show, this is not the case. 


When one solves (13), e.g., for negative curvature, /C = —r one obtains 


= “2 ) + >"’■• 


( 25 ) 


and evidently the associated and zero. 



T 


T 



X, Y 




X, Y 


FIG. 1: For Tifiat (figure on the left) the spacetime is essentially the same as Minkowski. 
Modulo an inessential singularity at the light cone, the causal structures are identical; 
hence, the spacetime is three-dimensional and flat all the way. In the case corresponding 
to Ticurved (figure on the right), we see that there is a change of signature crossing the 
light cone, that signals a behavior similar to a black hole horizon. In this latter case, the 


true spacetime is one dimension smaller, due to the constraint (27) that coordinates need 
to satisfy when curvature is present. 


Nonetheless, this time the metric ([^ reads 

r2 

9i.u{Q) = = 


J'2 _ _ y2 


1 0 0 
0-10 
0 0 -I 


(26) 


and the light-cone becomes essentially singular due to the change of signature, 
typical of horizons in black hole spacetimes. Hence, the light-cone is a proper 
Killing horizon, as shown already when considering Hawking effects reproduced on 
graphene; see, e.g., m- Let us now see what happens to the length of a position 
vector, 

\\Qf = 9,uQ^Q'' = r^>0. (27) 


That is a dramatic difference with the previous case (16) (amusingly, the difference 
comes about because of the different multiplicative constant factors: —1 for the flat 
vs —1/2 for the curved case). In the curved case the coordinates obey a constraint; 
hence, the effective theory lives in one dimension less, with the coordinates that we 
call P", a = 0,1. The behavior of the flat and curved case is depicted in Fig. 
Furthermore, the 6^ in the expression for the spin-connection (19) is not the dreibein 


anymore: the manifold index jj, now refers to a true curved space. All of this makes 
us conclude that the action that corresponds to the curved case (in the negative 
curvature case) is 


•^eff 


+ ^a)X, 


(28) 
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where everything has been dimensionally reduced to the two dimensions and, in 
particular, the metric is the induced curved metric: 

(2).px_ 

daisi ) dP°‘ ^ 

As well known, this two-dimensional theory has a quantum Weyl anomaly. Hence, 
interestingly, through the classical constraints we have a manifestation of a quantum 
Weyl anomaly. 

The Weyl (trace) anomaly is known to be in one-to-one correspondence with the 
Hawking radiation [13]. This is an interesting road to pursue for an alternative 
computation of the Hawking phenomenon on graphene as the one present in the 
literature mi- 


4. EQUIVALENCE OF THE STATIC HAMILTONIAN AND THE FULLY 
RELATIVISTIC APPROACHES 


To make contact with the literature focused on the phenomenology of graphene 
(see, e.g., [5]) we need hrst to move from the action in the coordinates to the 
Hamiltonian in the coordinates, keeping the curvature radius r finite. For E = 
— I ln(T^ — — Y^) -|- Inr = —t/r, the coordinates are 


T = 

X — cos x 

y = e'/V^Wsini 


(30) 

(31) 

(32) 


To write the Hamiltonian in the q coordinates, we have to consider the action (20), 
Legendre-transform it 


A 


dtdxdyC{q) 


jitiMy 


dQ 

dq 


{'K^dT'ilj){q) - H 


(33) 


and follow the steps presented in the Appendix A. This gives, for if = J dxdy'H{q), 


H = —ih / dxdy 'ip' (r* [vF{o'{y))]d d-, -f xpcp + vpPAi j xp 


where r* are the Pauli matrices, we have re-introduced h and vp, and 

d) = - 


(34) 


(35) 
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n3(J 


Al = — 


O’i'i 


n3(J 


r 


cosx , ^2 = — 


(J n, 


r + r2 


sinx 




^25 ^2;y 


with 


= 

= 

V2x = 

= 


r ey ay 

+ 5 


sin a; 


— a/ e^o- -f 7-2 cosx 


r e*" a,-, 


+' 

r 


cosx 


— a/ sin£ 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


The computations here were carried on for the S in (25), for which we can present 
the coordinates (30)-(32); hence, the expressions for 0, Ai and [vF]i ^ depend of that 


choice. Nonetheless, even though the detailed expression of those quantities change 
for S in (14), the structure of the Hamiltonian (34) remains the same for the case 


of interest of pure strain. 

We see here that, through this top-down method, we were able to reproduce all 
the terms except one that customarily appear in the literature of strained graphene; 
see, e.g., |5]. The latter is the one gauge field that gives unambiguous physical 
effects, and that couples to the spinors with an imaginary factor (an instance that, 
on its own right, is an indication that such field cannot be a Weyl held, see, e.g., 
ED- In the next Section we shall extensively comment on this held. Before moving 
to that, let us get one step closer to the language usually adopted in the graphene 
literature. 

In fact, the former expressions are written in the language of conformal factors a 
for the membrane and in isothermal coordinates whereas the usual approach employs 
strain tensors Uij = l/2{diUj + djUi) (where, as customary, the Ui measures the 
departure from the unstrained position), and Cartesian coordinates. 

Although we started oh with a fully relativistic formalism, due to the structure 
of the metric (|^, everything of the previous expressions necessarily depends only 
on the spatial coordinates. Henceforth, we can focus on the spatial metric only and 
make the customary ansatz that, in Cartesian coordinates {x,y), 


gij{x,y) ~ Sij + 2uij . 


( 42 ) 


On the other hand, this metric is related through a coordinate change to the one in 
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• = sublattice La o = sublattice Lb 


FIG. 2: The honeycomb lattice of graphene, and its two triangular sublattices, with our 
choice of the basis vectors, ai = |(\/3,3), 0,2 = |(\/3, —3), and si = £(0, —1), S 2 = 
|(\/3,1), S 3 = |(—\/3,1), with £ the carbon-to-carbon distance. 


0 




— T .. p2o-((}(g)) 


(43) 


This needs to be considered in its infinitesimal form, i.e., q^{q) ~ g* + m*, so that 
dq^/dq^ ~ <5* + dju\ which gives, at first order, Lij ~ 6ij + 2 -Ujj; hence. 


gij(x, y) ~ {5ij + 2uij) (1 + 2cr) . 


(44) 


Comparing the two expressions ( |42[ ) and (44), we obtain the wanted link between 
conformal factor, Cartesian strain tensor, and isothermal strain tensor 

Uij = a5ij + Uij . (45) 


5. THE HONEYCOMB STRUCTURE GAUGE FIELD 


Near the Dirac points (with our choice of the basis vectors in Fig. 2 = 

(±^,o)), we can describe the behavior of the tt electrons of graphene with the 
following Dirac Hamiltonian, 


H = —ivp 


cPx • V'0_ j , 


(46) 
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where V'lt are two-component Dirac spinors, a = ((Ti,ct 2 ), and a* = (di, —d 2 ), with 
dj the usual Pauli matrices, all emerging from the honeycomb lattice structure of 

Fig.|^ 

For strained graphene, the slope of the Fermi cone is changed and the Dirac 
points are shifted [H], + w±, where = ^Ai and w ±2 = so that 

the Hamiltonian is now 

H = -ivF j (fx ■ (^V + iA^ - ^pla* ■ (v - iA^ , (47) 

where Ai ~ {un — U 22 ) and A 2 ~ M 12 , or 

, (48) 


where z, j, k,... = 1, 2. This is the usual pseudogauge held interpretation of the strain 
effect in the graphene sheet [2]. The third rank anisotropic tensor contracted 
with the strain tensor Uij, is dehned in Appendix B, and will be analyzed in the 
following, as this will help to clarify the peculiar (anisotropic) nature of this gauge 
held. 

For homogeneous strain, the unstrained fermion pjo is related to the strained one 
'ip± via a phase transformation 

(49) 


The modihcation of the Hamiltonian (47) in the presence of an inhomogeneous strain 


is computed in [7]. In this case (49) becomes 




(50) 


where we follow the Dirac prescription [1SIE6], 'ipo is the unstrained solution. As a 
consequence of time-reversal symmetry, each Dirac point is charged with opposite 
sign. We can select a point xq as a reference strain, i.e., a point where the strain 
ehect w{xo) could be gauged away. The phase acquired in (50) could be interpreted 
as the circulation from the zero strain region xq to the strained one x. We stress on 
the fact that the solution (50) is formal: once the curl of w is non-zero, the integral 
in (50) is path-dependent (is a non-reversible process). Once we realize that the 
effect of the inhomogeneous strain is to shift the momenta of the Fourier modes in a 
space-dependent way, we can interpret the situation as the fermion is in presence 
of an effective gauge held A, which could give us a non-zero effective magnetic held 
if the curl of A is non-zero, giving rise to the characteristic Landau levels. This 
gauge held, frequently called in the literature “pseudo gauge held”, is the one we 
announced in the preceding part of this paper. Clearly, it could not have guessed 
from the QFT in curved space description. 
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y / ^ \ \ \ \ 

- ^ y / / 1 



FIG. 3: Strained graphene with a deformation vector defined as ft = {2xy,x‘^ — y^)uQ/AL. 


We can follow the trail of the non-trivial behaviour of A as the result of the 
contraction of the strain tensor Uij with the third rank tensor This tensor is 

nonisotropic^ and its presence is due to the triad {si, S 2 , S 3 } specific to the structure 
of graphene, which is built from two sub-lattices [2]. The nontriviality resulting 
from this contraction could be seen in a simple example [IS] . 

Consider the deformation vector u = {2xy, x‘^—y'^)uo/L, where mq is the maximum 
value of the strain and L is the length of the graphene sample, as represented in 
Fig-i This is a nonsingular vector held, in the sense that its curl and divergence 
is zero, so it is an irrotational vector held without any source or sink. However, we 
notice here that the vector held A resulting from the contraction of the corresponding 
strain tensor with jg ^ clockwise rotational vector held with constant curl, as we 
can see in Fig. This means that not only does the displacement vector u matter 
but also the orientation of this vector with respect to the near-neighbor vector triad 
{si}. So, the carries some memory of the lattice structure (short-distance 
behavior), even in the continuum limit (large-distance behavior). 

To see more concretely that this is the case, suppose that we have two unstrained 
graphene sheet samples. The honeycomb orientation of both samples are such that 
one of them, say it sample 1 , has the vector si of the triad parallel to the |/-axis, as in 
Fig-i while sample 2 has vector si = f'(—1, 0) parallel to the x-axis, see Fig. The 
orientation of the triad are both the same. A direct computation shows that the 
tensor are different for sample 1 and sample 2. Now let us apply to both samples the 


^ In fact, the only isotropic tensors of rank three are proportional to the Levi-Civita antisymmetric 
tensor eijk, which is zero in two dimensions [17j . 
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FIG. 4: Pseudogauge field produced using the deformation vector u = {2xy, — y‘^)uo/AL. 


same strain deformation vector shown in Fig. For the sample 1, the psendogange 
vector held gives us a constant pseudomagnetic held \B\ = jiUq, while for the 
sample 2 the pseudo-magnetic held is zero. So, even though in the continuum limit 
the two samples only diher by their orientation, hence should be indistinguishable 
by isotropy, the pseudo-magnetic held keeps track of the honeycomb orientation at 
the lattice level. A detailed study of the so-called “memory tensors” in the graphene 
honeycomb and Kagome lattices is found in 


Once the curl of w is not zero, equation (50) is not single-valued, as we can see if 


we take a loop around the origin. We are assuming that the variation of the strain 
tensor Uij is very small so, following the lines of [201121] , we can envisage a process 


to extract a physical meaning of the solution (50), even if it is not single-valued. 


Consider a small planar box on the graphene sheet (but very large compared with 
the lattice length t) situated at R (see Fig. |^. In the case of unstrained graphene 
{A = 0), the solutions for the 7r-electrons have the form^ '0j(r — R) = HlO >. Now, 
the idea is to strain the graphene sheet in such a way that in the box the strain is 
homogeneous (constant shift w of the Dirac point as in the previous section) and 
the associated magnetic held is almost zero in that region. The solution of the tt 


electrons in the box could be written using the Dirac prescription (50) 


K) = R)- 


(51) 


Because we are in the second quantization formalism, tp are operators and the wave packets are 
these operators applied to the vacuum. 
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FIG. 5: Two samples of graphene unstrained sheets. Sample 2 is clockwise rotated respect 
to sample 1 by 7r/2 radians. Once we apply a deformation vector shown in Fig. the 
sample 1 acquires a constant pseudo-magnetic field |.B| = ^uq, while for the sample 2 the 
pseudo-magnetic field is zero. 


We can see that "0 in (51) is single valned in r and R locally. Now let the box 


be transported aronnd a loop C. After completion of the loop, there will be a 
geometrical phase change that can be compnted nsing Berry’s formnla 


t(C) = ! (t R)Vgiir- R) ■ dR 


(52) 


C 


Taking into acconnt that the "0 are normalized, we end np with the resnlt 


7(C) = - J MR) ■ dR = -4, 


(53) 


where is the magnetic flnx enclosed by the loop. 


6. CONCLUSIONS 

When only strain is present, the only gange field that has nnambignons physical 
effects, is the one jnst discussed in Section The structure we saw there, is rem¬ 
iniscent of the gauge field arising in the Aharonov-Bohm effect, although strictly 
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FIG. 6: Loop C of the box (in blue) situated at R. All the fermions inside the box could 
be described as R). 


speaking, in order to have this effect, we need a magnetic flux line crossing the loop 
C and a zero magnetic field outside the line. This is also reflected in the fact that the 
flux obtained in (53) is not constant and depends on the area enclosed by the loop. 
However, it could be possible to imagine a strain vector concentrated in some region 
of the graphene and assumed to be zero outside this region. Then, this procedure, 
recover the Aharonov-Bolim result. The Aharonov-Bohm effect is not new in the 
graphene literature and some examples of strain vectors were proposed in order to 
see this effect with particular procedures, see for instance |23l 1^ . 

This gauge field could not have guessed from a top-down approach based on the 
standard QFT in curved spacetime. In fact, we saw here that if we define standard 
classical functions on graphene membrane the only fields emerging are the metric 
and the spin connection, which do not reproduce the behavior of this pseudogauge 
held. This U{1) held can be put in correspondence with quantum held theoretical 
structures, such as the quantum anomalies. To see it, one just needs to realize that 
the origin of such gauge held is entirely quantum mechanical, and related to the fact 
that a “translation” in conhguration space T ■. x ^ x + u, that is the straining of the 
graphene membrane, is necessarily associated with a “translation”® in momentum 
space B : k ^ k + v. Hence, in the hrst quantization language of the wave functions, 
those operations are carried on by quantum operators W : T x H —)■ C that, for the 
very meaning of quantum mechanics (that is the Heisenberg uncertainty principle) 


® A better name for this operation is ’’Galilean boost”, that is why we use “B” for it. On this 
note see, e.g., pS] . 
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need to obey 


W(ri, 5 i)W(T2, B 2 ) = 52)W(Ti, Si) , (54) 

that are recognized, in the second qnantization language, as instances of the nontriv- 
iality of the quantum field theoretical vacuum, and in turn of the quantum anomaly 

izg. 

Therefore, our simple table-top laboratory, can help explore yet another arena of 
fundamental physics, that is the deep meaning of the quantum anomaly®. 

On the other hand, this pseudogauge held, for the use of graphene we have in 
mind, is intriguing also because, as explained in the paper, it is a memory of the 
lattice structure (that is the physics of the wavelengths comparable to £, the lattice 
spacing) in the continuum limit (that is the physics of the large wavelengths). That 
means that it is a relic at low energy of the high energy behavior of the system. This 
is yet another reason why we cannot reproduce this held from a top-down approach, 
because in such approach the isotropy of the graphene membrane is tacitly assumed. 

Effects of this kind would be of paramount importance to bring high energy 
theoretical constructions under the control of experiments. One example that comes 
to mind is the standard model extension of |28] , where tensorial helds that are relics 
of the Lorentz invariant high energy string theory combine with the helds of the 
Standard Model (SM) and their derivatives within Lorentz violating terms that 
have the form 

^(SM helds and derivatives)^"'^ , (55) 

where ~ ^pianck- “graphene universe” ^pianck ~ 

We plan to explore these scenarios, and to address the full variety of possible 
gauge helds in graphene, namely those arising beyond the pure strain limitation, in 
later work. 
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® In our view, the very existence of a quantum anomaly is a sign of lack of understanding of how 
nature works at the most fundamental level, or, in other words, it is a sign of our ignorance of 
what is a proper quantum Noether theorem [27] . 


18 



Appendix A: Prom (2 +l)-dimensional Lagrangian to static two-dimensional 
Hamiltonian 


To be specific, let us focus on the case on negative constant Gaussian curvature, 
although similar formulae must hold the pure strain/fully flat case. The coordinates 
are 


T = +r2 

A = e*/'’e'"(^)cosx 
Y = sin x 


(Al) 

(A2) 

(A3) 


and 


S = -- \n{T^ - - Y^) + \nr = -t/r 


To write the Hamiltonian in the q coordinates, we have to consider the action (20) 

A = i f dTdXdYe^^ + Hy 


We have to write this in the coordinates. We need the Jacobian 




dq'' 


•y/g2(T(y) _j_ J..2 


(A4) 

(A5) 


the three terms 


YjJ' — 

Sv = 

Ey = 


T 2 - A 2 - y 2 r 

X e-*A 


J'2 _ Ji^2 _ y2 ^2 

Y _ e-*A 

_ J ^2 _ y 2 ~ J.2 * 


t/r - 

— V e2°'(5) -|- r2 

(A6) 

COSX 

(A7) 

sin X 

(AS) 


and to re-express the derivatives, e.g., dx = Xxdx + Vxdy, etc., where, as 
(Tyijj) = dya{y), xx = dx/dX, etc. Then we use 


A= dtdxdyC{q) = / dtdxdy 


dQ 


dq 


{Ti^dT'ilj){q) - U 


usual, 

(A9) 
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from which we can read off the Hamiltonian H = J dxdy'H{q). The final expression 

is 


H = —ihvp / dxdy 


r ay 


+ 


- 0 ' 


\/e 2 ^ +' 


V’ 


(AlO) 


+ 


^ — e sin X + 


e ^ 


ay 


cosx dy ^^ cosxj tp (All) 


+ (e cos X dx + 


e 6^°" + 


cTy 


sin X dy -\—^ sin x) 'ip 


(A12) 


where we have reintrodnced h and the Fermi velocity vp- Note that nothing depends 
on t, as it must be. 


The formula above for H gives a field of the type 'ip'^'ip in (AlO) 

1 


= - e^'^o-y , 
r 


(A13) 


the Ai and A 2 fields in the non-derivative terms of (All) and ( A12[ ) (here we identify 
Ai and A 2 according the Pauli matrices and r^) 


n3(J 


n3(J 


Ai = — 


COSX , A2 = — 


r y'g2cr ^2 ’ r -^4(2^^772 

and the space-dependent, inhomogeneous Fermi velocity tensor 


smx 


(A14) 


[■<'f(ct(»))1, ' = vf 


Vlx Vly 
V2x V2y 


with 


Vlx = 


r e" ay 

+ '1 


smx 


Vly 

V2x = 


=- \/cosx 


r e" a,-. 




cosx 


V2y 


g'-' 

y = -\/-1- sinx 

r 


Therefore, obtaining the expression (34) we have in the main text. 


(A15) 

(A16) 

(A17) 

(A18) 

(A19) 
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Appendix B: Some details of the tight-binding computations 


We start with the effective function of the hopping energy respect to intercarbon 
distance na, 


ti ^ -?;exp 


- 1 ) 


(Bl) 


where rj ~ 2.8 eV is the equilibrium hopping energy and s', with i = 1,2, 3, are the 
variation of the basis vectors s^. In order to deal with inhomogeneous strain, we 
consider the expansion of (Bl) up to the first derivative of the strain tensor Uij{x), 
i.e.. 


ti = r] 


1 - 




(B2) 


where with some abuse of notation, (sj)™' stands for the m component of the vector 
Si, the indices k,m,n are contracted (dummy indices) and dk = gf*. The tight- 
binding Hamiltonian, in the second quantization formalism^, could be written as 


i=3 

H = — (^a^{x)tib{x + Si) -1- c.c.) , 

x&La *=1 

where La stands for the sublattice A (see Fig. and, in order to have a Hermitian 
Hamiltonian, we added the complex conjugate of the first term. Using the following 
expressions for the Fourier transforms 

a(£) = ^e“aj, • «’"(£) = • 

k k k 


we can write the Hamiltonian in the Fourier space as 


H = 



0 


i,q 


i,q 



(B3) 


where we used the symmetry property of Uij and defined as 


= 5(g) + ^(si)”'Mm,g(si)” - (si)^(si)”"(si)”gjgmM 


n,q' 


^ If we use the second quantization formalism, some subtitles appear with respect to the vacuum 
where the creation and annihilation operators act. This subtitles will not be considerer here 
because the system is simple enough but, in the case of curved graphene, extra care must be 
taken. This is because the presence of defects on the honeycomb makes the vacuum nontrivial 

m- 
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Now, we expand the Hamiltonian aronnd one Dirac point K±. For the sake of 
simplicity, we expand aronnd Kj^ = 0 j, snch that k = K+ + p. We can work 


a little bit more in the matrix content of (B3) as 


n T'^ \ j / n 

- *’'? I = 7^ ■ '^<^3 (I + — <?)) I A’^ m 

0 J i \ 0 

where in the last eqnality we made use of the identity® [7] 

0 e-hF+-<jO-si 




0 


- ' ^2^3 • 


At this point, it will be useful to show the following identities [7j 


i=3 




2=1 


1 Q 

^ \ ^^ \m(^ \n ^ ^mn 

2=1 

2=1 

1 3 

• 1 ^ 

2=1 

The tensor defined in the third identity, is an invariant under the discrete C 3 

rotations. This tensor is very important in the discussion about the anisotropy of 
the strain-induced gauge field. For our choice of basis vectors {^}, its only nonzero 
components are 7^222 _ _/^i 22 _ _/^ 2 i 2 _ _/^ 22 i _ other 

three tensors are all isotropic. Making use of all this, doing some standard algebra 


In the case of expanding around the other inequivalent Dirac point 


K_ 


formula is 


C O 

gi(iA-(3)-sl Q 
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and going back to the confignrations space via the anti-Fonrier transforms of a, b 
and M™', we end np with the following Hamiltonian for the inhomogeneons strain 

H = {vjyndm + ivpA^ - ^fF^) -0+, 

X 

where Vjm = vp (^5jm — ^{unnbjm +is the celebrated space-dependent 
Fermi velocity [7], is the psendogange held, and F-^ = 

I {dmUjm + \djUmrn) IS a connection-like coefficient. The Fermi velocity in nn- 
strained graphene is vp = ^rji. 
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